Matter power spectrum for the generalized Chaplygin gas 
model: The relativistic case 
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Abstract 



The generalized Chaplygin gas (GCG) model is the prototype of a unified model of dark 
energy (DE) and dark matter (DM), ft is characterized by equation-of-state (EoS) parameters 
A and a. We use a statistical analysis of the 2dFGRS data to constrain these parameters, 
tn particular, we find that very small (close to zero) and very large values (a 3> f ) of the 
equation-of-state parameter a are preferred. To test the validity of this type of unification 
of the dark sector we admit the existence of a separate DM component in addition to the 
Chaplygin gas and calculate the probability distribution for the fractional contributions of 
both components to the total energy density. This analysis favors a model for which the 
Universe is nearly entirely made up of the separate DM component with an almost negligible 
Chaplygin gas part. This confirms the results of a previous Newtonian analysis. 
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The Chaplygin gas model represents the best known proposal of a unification of DM and DE 
into a single fluid pQ . In its generalized form [2] , the Chaplygin gas is a fluid with an equation of 
state p c = —A/pc, where p c is the pressure and p c is the energy density. Solving the conservation 
equation p c + 3a/a (p c + p c ) = 0, where a is the scale factor of the spatially flat Robertson- Walker 
^ ■ metric, the result for the energy density is 



Pc 



A + Ba~^ l+C 



V(i+") 

(1) 



The present value of the scale factor is normalized to do = 1- Combined with Friedmann's 
equation and assuming a > — 1, the solution ([I]) interpolates between a matter phase for a <C 1 
and a de Sitter phase for a>l. In this sense, the GCG can play the roles of dark matter in the 
past and of dark energy at present and future times, thus unifying the two components of the dark 
sector of the cosmic substratum into a single component. The original Chaplygin gas corresponds 
to a = 1. It can be traced back to the Nambu-Goto action in the light cone parametrization [3], 
that is, it can be seen as rooted in a string theoretical context. The generalization to a ^ 1 is 
phenomenological [2j. However, the GCG model has the advantage of providing an interpolation 
between the ACDM model, represented by B = and a = 0, and other models, like, e.g., viscous 
models [HOSE]- Similar approaches are the so-called Cardassian models 0|8]. Thus, the GCG 
model covers a large variety of dark energy models. 
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A great deal of effort has been made in order to constrain the different parameters of the GCG 
model, mainly the parameter a. This includes data from supernova type la [9], the anisotropy of 
the cosmic microwave background |10j . baryonic acoustic oscillations the integrated Sachs- 
Wolfe effect |12| and the matter power spectrum (see |13j and references therein). Constraints 
from combined data sources have been obtained in |14| and [15] , 

The supernova type la data were shown to be well described by the GCG gas, which seems 
to be favorable for the unified model. But the situation changes if the model is confronted 
with matter power spectrum data. In this brief communication we generalize our previous 
Newtonian analysis |13| of the 2dFGRS program [16J to the fully relativistic case. We consider 
three configurations of a Chaplygin-gas cosmology 

(i) At first we assume the material content of the Universe to consist of a GCG, radiation and 
a pressureless baryon component with a fraction of Vt m Q = fJ&o = 0.043, as suggested by the 
five-year WMAP data. In other words, we prescribe a unified model of the dark sector. The 
result is that, together with very small values of a, for which the GCG behavior is close to that of 
a cosmological constant, also values a > 2 have a high probability. This confirms both an earlier 
qualitative but gauge invariant analysis of matter power spectrum data [17J and our previous 
simplified Newtonian analysis [13J. Still more surprising is that very large values, i.e. values of 
several hundreds for a, are preferred by the data. This confirms an independent investigation of 
the integrated Sachs- Wolfe effect with preferred values up to a ~ 350 [18J. 

(ii) At second we consider the original Chaplygin gas, i.e., the case a = 1, but we leave the matter 
fraction free, thus admitting that the matter component is not just made up by the baryons. This 
is equivalent to allow for a separate DM component in addition to the GCG. For this reason we 
write the fractional density parameter for the pressureless matter today as f2 m o = 0.043 + QdmO- 
This additional freedom is used to test the unified model of the dark sector. The unified model 
can be regarded as favored by the data if the probability distribution function (PDF) for the 
matter fraction is large around the value that characterizes the baryon fraction. If, on the other 
hand, the PDF is largest at a substantially higher value, the unified model has to be regarded 
as disfavored. Our fully relativistic study confirms the results of a previous Newtonian analysis 
according to which the matter fraction does not only not peak around the value for baryons, 
but it is maximal at values of f2 m o close to f2 m o = 1- in other words, the data seem to prefer 
a model which is close to the Einstein-de Sitter universe. This contradicts the results from the 
SN la data according to which the unification scenario, i.e, a universe filled almost exclusively 
by the Chaplygin gas component, is favored [9]. 

(iii) At third we extend the analysis of (ii) to a / 1. As far as the value of f2 m o is concerned, we 
recover the result of (ii). The behavior of a is similar to that described under (i), although here 
it is less important since the GCG is subdominant. 

Strictly speaking, models with a > 1 seem to be un-physical since they correspond to a 
superluminal sound speed. However, with certain modification causality can be preserved |17] , 
To see the influence of a large a on the background dynamics of a GCG cosmology, we consider 
the EoS parameter 

Pc A A 

— = — = ? ^ 57i — 7 with A = — -T . (2) 

Pc A+(l-A) a- 3 ^) p^ 1 { ' 

The parameter a influences the value of a at which the transition from decelerated to accelerated 
expansion occurs. Assuming a unified model and denoting the transition value of the scale factor 
by a q , one has 

Pc 1 A 



Pc~ 3" A+{l-A)a- 3il+a) ' 



(3) 



Solving for a q yields 



IT ' (4) 



Since < 1 this means, that for growing values of a the transition period a q approaches 

a = 1, i.e., the present epoch. For very large a, the matter period is longer and the transition 
to accelerated expansion occurs suddenly and more recently than in the ACDM model. This 
property has also been discussed in |18j . 



2 Basic set of equations 

Our starting point are Einstein's equations coupled to a pressureless fluid, radiation and to the 
GCG. They read, 

rpfll/ _ Q rpflU _ Q rpflV _ q 

The superscripts (subscripts) m, r and c stand for "matter", "radiation" and "Chaplygin". We 
assume a perfect fluid structure for the cosmic medium as a whole and also for each of the 
components, 

T>» = pu»u v + p {sT ~ u»u v ) , Tjf = p A u\u\ + p A {gT ~ u\u\) , A = m, c, r . (5) 
Using now the flat Friedmann-Robertson- Walker metric 

ds 2 = dt 2 - a(t) 2 [dx 2 + dy 2 + dz 2 } , 
and identifying all the background 4-velocities, Einstein's equations reduce to 

d\ 2 8vrG 8ttG 8ttG 

-Pm H ^— Pr H 5— Pc, (6) 



a J 3 ' 3 r 3 

2-+f-V = -8vrGp c , (7) 

Pm + 3- =0 => p m = p m0 a , (8) 

Pm 

Pr + 4- = => p r = PrOO" 4 , (9) 

d r b iV(i+«) 

p c + 3-(p c + p c ) = (p e = -A/p?) => Pc =U + -^— j \ . (10) 

The perturbed equations in the synchronous gauge take the form, 

- + -h-47rG(5p + 3<5p) = (11) 
2 a 

6p + ^(5p + Sp) + (p + p)(e-^\ = 0, (12) 
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where p and p stand for the total matter and pressure, respectively, and 9 = 5u l i . 
In terms of the components, we end up with the following equations: 

^ + %-ATTG[5p m + 5p c + 5pr + 2>{5p m + 5p c + 5p r )] = 0, (14) 

5p m + ^ {5Pm + 5Pm) + (Pm + Pm) ( m - ~ J = 0, (15) 
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with 9 m = 5u l mi , 6 C = 5u l ci and 9 r = 5u l ri . 
With the definitions 
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the set of first-order equations becomes 
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where 
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and feg 1 = 3000 /iMpc. 



3 Bayesian analysis 

The matter power spectrum is defined by 



where 5k is the Fourier transform of the dimensionless density contrast 5. We will constrain the 
free parameters using the quantity 

where V° is the observational value for the power spectrum, V\ is the corresponding theoretical 
result and <7, denotes the error bar. The index i refers to a measurement corresponding to 
given wavenumber. The quantity (|34p qualifies the fitting of the observational data for a given 
theoretical model with specific values of the free parameters. Hence, x 2 is a function of the free 
parameters of the model. The probability distribution function is then defined as 

F(x n ) = F e-x 2 ^ 2 , (35) 

where the x n denote the ensemble of free parameters and Fq is a normalization constant. In 
order to obtain an estimation for a given parameter one has to integrate (marginalize) over all 
the other ones. For a more detailed description of this statistical analysis see Ref. [9]. In the 
evaluation of the probability distribution, we will use the total x 2 - With the prior of a spatially 
flat universe, the analysis has to take into account three free parameters: a, f2 m o and A (or 
equivalently Qdmo) an d A). As already mentioned, we will consider three cases. For case (i), 
the unified model, the fraction £l m Q is fixed to m o = £lbo = 0.043 and there are only two free 
parameters, a and A. The PDFs are visualized in figure [TJ The upper left and central panels 
show the two-dimensional PDF for A and a for different ranges of a. The lighter the color the 
higher the probability. In the upper right panel the one-dimensional PDF for A is depicted. It 
has maxima close to zero and close to unity. Recall that A = means a pressureless medium 
while A = 1 corresponds to a cosmological constant for any value of a. Both lower panels show 
the one-dimensional PDFs for different ranges of a. There is a maximum at a = which, in the 
background, corresponds to the ACDM model. Also values a > 2 have a high probability. This 
feature is both in agreement with a previous gauge-invariant perturbation analysis [17J and with 
our simplified Newtonian model |13] , The more surprising result, however, is the existence of 
a maximum at very large values of a, as can be seen in the lower right panel of figure [TJ This 
backs up an independent study of the integrated Sachs- Wolfe effect |18| with values of a up to 
350. Here, we obtain a maximum for a = 240. 

Case (ii) is the original Chaplygin gas with a = 1 and we have two free parameters as well, 
O m o and A. The results of the statistical analysis are visualized in figure [2j In the left panel 
the two-dimensional PDF for Q m o and A and a is seen. The lighter the color the higher the 
probability. The center panel shows the one-dimensional PDF for A. It is maximal near A = 1. 
The right panel depicts the one-dimensional PDF for Sl m o- It is maximal around £l m o ~ 0.95, a 
value close to unity, i.e., much larger than the baryon fraction. According to our criterion, the 
unified model is clearly disfavored, since the GCG fraction is of the order of 5%. 

In case (iii) all the three parameters a, £l m o and A are left free. The two-dimensional PDFs 
for all binary combinations are shown in figure [3l the one-dimensional PDFs in figure [H Again, 
the maximum value of Cl m o is close to unity, confirming the result of case (ii) according to which 
the contribution of the GCG component is almost negligible. The behavior of a is similar to 
that of the unified model of case (i), although this result is less important here since the GCG 
is subdominant. 

We remark that a more detailed, quantitative estimation for the parameters of the model 
is doubtful due to the existence of different regions of high probability implying that the one- 
dimensional PDF is not gaussian. 

4 Conclusions 

Our fully relativistic analysis corroborates the results that were previously obtained within a 
simplified Newtonian analysis [13J. As in the Newtonian case, the PDF of A is sensitive to the 
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Figure 1: PDFs for the unification scenario (f2 m o = ttbo = 0.043). 




Figure 2: PDFs for the original Chaplygin gas case a — 1). The point in the figure in the left indicates the 
maximum in the two-dimensional distribution. 

number of free parameters considered. For a there are two regions of high probability: one region 
near zero and a second one for a ~ 240. Although the existence of this second region comes as 
a surprise and represents a new result compared with |13j . it confirms an entirely independent 
investigation, using the integrated Sachs- Wolfe effect [18J. The probability distribution for fl m o 
is essentially the same as in the Newtonian case: a universe with negligible Chaplygin gas 
component is favored. In fact, our analysis of the large-scale structure data from the 2dFGRS, 
taken separately, does not favor an accelerating universe. This is in striking contrast to the 
results for the homogeneous and isotropic background on the basis of the supernova data [9]. 
This may imply that the Chaplygin gas model should be discarded. But on the other hand 
it is well known (see also |13j), that the matter power spectrum generally does not sufficiently 
constrain the dark energy component. Further studies, in particular the crossing with other, 
independent tests are desirable to assess the status of Chaplygin- gas cosmologies. 




Figure 3: Two-dimensional probability distribution for the different combinations of the parameters a, SldmO 
and A. 
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Figure 4: One-dimensional probability distribution for a, QdmO and A. 
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